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Abstract 


We investigate the chemical transition of simple molecules like C2H2 and HCN into 
aerosol particles in the context of Titan’s atmosphere. Experiments that synthesize 
analogs (tholins) for these aerosols can help understand and constrain these polymer- 
ization mechanisms. Using information available from these experiments, we suggest 
chemical pathways that can link simple molecules to macromolecules, that will be the 
precursors to aerosol particles: polymers of acetylene and cyanoacetylene, polycyclic 
aromatics (PAHs), polymers of HCN and other nitriles, and polyynes. Although our 
goal here is not to build a detailed kinetic model for this transition, we propose param- 
eterizations to estimate the production rates of these macromolecules, their C/N and 
C/H ratios, and the loss of parent molecules (C2H2, HCN, HC3N and other nitriles, 
CqH 0) from the gas phase to the haze. We use a 1 -dimensional photochemical model 
of Titan’s atmosphere to estimate the formation rate of precursors macromolecules. 
We find a production zone slightly lower than 200 km altitude with a total produc- 
tion rate of 4xl0~ 14 g cm -2 s _1 and a C/N ~ 4. These results are compared with 
experimental data, and to microphysical models requirements. The Cassini/Huygens 
mission will bring a detailed picture of the haze distribution and properties, that will 
be a great challenge for our understanding of those chemical processes. 
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1 Introduction 

The dense atmosphere of Titan is composed of approximately 98% nitrogen molecules and 2% 
methane. Photochemistry occuring in this atmosphere yields a suite of other gaseous compounds, 
including hydrocarbons and nitriles (e.g. ethane, hydrogen cyanide). An aerosol layer, which hides 
Titan’s surface in visible light, is also present in the stratosphere. The composition of the gas 
and the microphysics of the haze have been studied through observations (both ground-based and 
from the Voyager spacecrafts) and modeling. But the formation of this haze, from photochemical 
processes in the gas phase to the smallest aerosol particles considered in the microphysical models 
(which we will call “precursors”), is not well understood. Photochemical models (Yung et al, 1984; 
Toublanc et al., 1995; Lara et al, 1996; Lebonnois et al, 2001) do not follow the fate of gaseous 
compounds with more than 6 heavy atoms (C and N), while microphysical models (McKay et al, 
1989; Toon et al, 1992; Rannou et al, 1995, 1997) assume the existence of small particles (of radius 
around 13A) as the source of material, then follow the microphysical evolution of the particles while 
conserving mass. Only one work has been done to detail the chemistry of hydrocarbons up to Cgo 
macromolecules (Dimitrov and Bar-Nun, 1997), but this work does not include N-bearing species, 
and has not yet been applied to simulate Titan’s atmospheric composition. In Dimitrov and Bar- 
Nun (1999), the authors present a model of the agglomeration of aerosols that uses this study. 
Their model can well reproduce the laboratory synthesis of aerosol particles from Bar-Nun et al 
(1988). It was also applied to describe a possible scenario for the synthesis of Titan’s aerosols, 
but no direct comparison was made with observations or with other detailed microphysical models 
(e.g. Toon et al, 1992; Rannou et al., 1995, 1997). 

In this paper, we address the question of how to link photochemical models to microphysical 
models in a simple, parameterized way, without including a detailed kinetic model to describe 
all reactions involved in the macromolecules formation. In order to couple these two kind of 
models, it is necessary to quantify the production rate of precursors from photochemical reactions, 
and to evaluate the chemical loss of gas phase molecules to the haze. This does not concern 
the possible condensation of gas phase molecules on aerosol particles in the low stratosphere and 
troposphere (typically 40 to 150 km altitude). This is a problem that needs to be treated in the 
microphysical models, and that does not affect the source function of the haze. In the microphysical 
models, no addition of mass from the gas phase is considered, which means that chemical growth of 
macromolecules with addition of gas phase small molecules (e.g. C2H2, HCN, HC3N) is assumed 
negligible when the smallest aerosol particles are considered. We w r ill call this transition the 
precursor level. 

The principles adopted in this paper are the following: 
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• given a set of chemical species included in the photochemical model, the production of aerosols 
occurs through a set of initial reactions that will produce the first of the heavy molecules 
(not included in the photochemistry). These molecules will irreversibly yield aerosol material. 
Meanwhile, the growth of these macromolecules is a sink for gas phase molecules included in 
the photochemistry (e.g. C2H2, HCN, HC3N). 

• during the initial steps of aerosol formation, the macromolecules grow from the first heavy 
molecules by including small gas phase molecules. Nucleation eventually occurs. Particles 
grow in mass and radius through several processes: chemical additions (chemical activity 
continues after nucleation), transition of gas phase macromolecules to the condensed phase, 
coagulation. By the time the particles reach a radius of approximately 1-2 nm, microphysical 
considerations suggest that the further coagulation of particles occurs at constant mass, i.e. 
the addition of mass from the small gas phase molecules has become negligible. We define 
the precursor level as the transition between the regime where chemical addition of small 
molecules is significant and the regime where microphysical processes conserving the overall 
mass of aerosol material dominate. The source function used by microphysical models is the 
mass flux of aerosol material through this transition. 

• this aerosol material, at the precursor level, is composed of macromolecules that have a 
certain mean length, i.e. that have incorporated, in average, a certain number of gas phase 
molecules during their chemical growth. 

Based on these principles, we estimate the source function of precursors and their bulk elemental 
composition through these steps, using three sets of parameters: 

• we propose three pathways for the initial reactions and the growth of macromolecules (poly- 
mers of C2H2 and HC3N, PAHs that may include HCN and HC3N, and polymers of HCN and 
nitriles). These propositions are based on experimental results published by several teams 
(see Table I). 

• the reaction rate coefficients (k, in cm 3 s -1 ) for these initial reactions are a first set of 
parameters. 

• for each pathway, we propose gas phase molecules that may be incorporated in these macro- 
molecules. The proportions of each molecule in the global macromolecule are controlled by 
a second set of parameters. 

• the yield limiting the growth of a macromolecule (precursor level) is parameterized by the 
mean number of gas phase molecules incorporated in one macromolecule. We also refer to 
this parameter as the ‘length’ of the polymers. 

• the mass flux of aerosol material at the precursor level (i.e. the source function of the 
aerosols) is equal to the product between the mean molecular mass of macromolecules and 
the production rates of the first heavy molecules (initial reactions). 

Laboratory experiments have synthesized tentative analogs to Titan’s aerosols, in order to obtain 
their properties. The study of these products have yielded some hints about the structure of the 
constituent macromolecules of the tholin particles, and in some cases the C/N and C/H ratios 
(Table I). How close tholins represent Titan’s aerosols is in debate, but still, understanding the 
process which makes tholins can certainly help. Using these experiments, we have identified three 
possible pathways leading from parent molecules (C2H2, HCN, HC3N and other nitriles, CgHg) to 
macromolecules that will be the building blocks of the precursors. These pathways are described 
in Section 2. In Section 3, the parameterization of these pathways is detailed, and quantities of 
interest (mass production rates of the different channels, loss of the parent molecules, C/N and 
C/H ratios) can be evaluated. This model is then applied to the atmosphere of Titan in Section 4. 


2 Information from experiments 

Numerous experiments have been done in order to synthesize tholins, analogs for Titan’s aerosols. 
The conditions in the laboratories have been varied while trying to get as near to conditions in 
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Titan’s atmosphere as possible. Table I presents a summary of these experiments, with some 
results on the products’ analysis. Since the optical properties of tholins provide good matches to 
the properties of Titan’s haze (Khare et al., 1984; Sagan et al ., 1992), it is reasonable to assume that 
the production processes for tholins are close to those occuring in Titan’s atmosphere. Most of the 
C/N ratios measured in the tholins are low (between 1.5 and 11), so we expect a fairly important 
amount of nitrogen to be incorporated in the macromolecules during the chemical processes. 


Table I: Experiments of laboratory synthesis of Titan’s aerosols’ analogs 


Papers 

ES^ 

T(K) 

P(mb) 

Initial gas mix 


Aerosol particles 

C/N 

C/H 

Khare at. al., 1984 

D 

300 

0.2 

90% N 2 , 10% CH 4 

- 

tholins 



Sagan at al, 1984 

D 

300 

>10 

90% N 2 , 10% CH 4 

- 

tholins 

2 

0.6 

Bar-Nun at al, 1988 

UV 

300 

270 

Ar, 5% C 2 H 2 

- 

polyC 2 H 2 







Ar, 5% C 2 H 4 

- 

polyC 2 H 4 







Ar, 5% HCN 

- 

polyHCN 



Thompson at al., 1991 

CP 

300 

17 

90% N 2 , 10% CH 4 

X 

tholins 






0.24 






Scatterwood at. al, 1992 

UV 

300 

73 

N 2 , 2% C 2 H 2 

- 

polyC 2 H 2 







N 2 , 2% C 2 H 4 

- 

polyC 2 H 4 







He, 2% HCN 

- 

polyHCN 







Mix 

- 

polyC 2 H 2 



McDonald at al., 1994 

CP 

300 

1 

90% N 2 , 10% CH 4 

- 

tholins 

1.5 

0.6 

Coll at al, 1995 

D 

150 

900 

90% N 2 , 10% CH 4 

X 

tholins 

11 

1 

McKay, 1996 

D 

300 

1000 

90% N 2 , 10% CH 4 

- 

tholins 

5.5 

1 

Clarke and Ferris, 1997 

UV 

300 

330 

90% C 2 H 2 , 10% HC 3 N 

- 

polyC 2 H 2 /HC 3 N 

6.4 

1.6 

Coll at al, 1999 

CP 

150 

2 

98% N 2 , 2% CH 4 

X 

tholins (polyHCN ?) 

2.8 

0.8 

Khare at al., 2001 

D 

300 

71 

90% N 2 , 10% CH 4 

X 

tholins 




(a) Energy Source: D = discharge, CP = cold plasma, UV = ultraviolet light 


(b) An “X” in this column indicates that the composition of the resulting gas phase has been studied. 


We used the information provided by these experiments to identify several pathways from gas 
molecules to aerosol particles. Possible pathways that we will analyse in this paper include different 
kinds of polymers (C2H2, HC3N, HCN), polyynes and cyanopolyynes and aromatics (PAHs). 

Polymers of C2H2 and HC3N 

These are obtained by polymerization of acetylene, or of a mixture of acetylene and cyanoacetylene. 
One possible structure for these polymers is described in Clarke and Ferris, 1997. Clarke and Ferris 
also indicate the incorporation of other molecules in these polymers (CH4, C2H6, CO), but it seems 
to be of second order. This incorporation should also occur in the case of Titan’s atmosphere, but 
we will restrict our study to acetylene and cyanoacetylene only as a first approximation. 
Polyaromatics 

Polyaromatic compounds have been detected in laboratory tholins from the Khare et al. (1984) 
experiment (Sagan et al., 1993; Khare et al., 2001). Previous theoretical work has investigated 
the growth of the PAHs from benzene by addition of acetylene (Wang and Frenklach, 1994; Wong 
et al., 2000; Bauschlicher and Ricca, 2000). Nitrogen atoms could be incorporated in the cycles 
through the addition of HCN or HC 3 N instead of acetylene (Ricca et al., 2001). Reaction rates 
and pathways for this growth are mostly studied at high temperatures (combustion), but some 
studies cover as low as room temperature, which can allow a tentative extrapolation for Titan’s 
atmospheric conditions (100-200K). 

Polymers of HCN 

Coll et al. (1999) report that polymers of HCN could be part of the tholins they produced. In 
this experiment, as well as in other similar experiments (Thompson et al., 1991; McDonald et al., 
1994), nitrogen has to be incorporated in large quantities to get the observed low C/N ratio, and 
the obvious candidate is HCN, considering that it is the dominant N-bearing compound (other 
than N 2 ) in the gas phase. The question is, in these experiments, what exactly is produced? HCN 
could be incorporated as a copolymer in an aliphatic polymer, together with C 2 H 2 and HC 3 N (but 
we are unaware of any possible structure), or it could polymerize on its own. Matthews (1992) 
presents a study of these polymers, and Minard et al. (1998) discuss possible structures for their 
formation. Thompson and Sagan (1989) propose another possible structure for a more general 
polymerization of nitriles. 

Polyynes and cyanopolyynes 

The structure of these molecules is A-(C=C)„-B, with A and B being either H or CN. Many 
members of this group have been detected in the gaseous products (Coll et al., 1995; 1999). Only 
C 2 H 2 , HC 3 N, and C 4 H 2 have been observed in Titan’s atmosphere. C 4 N 2 has also been detected 
as a solid in the polar atmosphere of Titan. The C/N and C/H ratios in these compounds are 
high: the number of C atoms grows with the length of the molecule, while the number of H 
and/or N atoms is only 1 or 2. Their fate in Titan’s atmosphere can be photolysis to smaller 
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compounds, condensation in the low stratosphere, or incorporation in the aerosols. Since the C/N 
ratios available in the experiments are rather low (the largest is 11, in Coll ct al., 1995) and the 
C/H ratios are all around 1, the polyynes and cyanopolyynes formation and incorporation may 
not be a dominant growth mechanism for experimental tholins. Despite the difference between the 
pressures used in experiments and the pressures in the atmosphere of Titan (in the [200-600] km 
range in altitude), we assume here that these products are also of second order in the constitutive 
material of aerosol precursors. 

Whether the products obtained in the experiments are long chains or aromatics, or a mix of both, 
is in debate. In this paper, we will concentrate on the first three mechanisms: polymerization of a 
mixture of acetylene and cyanoacetylene (polymer 1); PAHs’ growth from addition of acetylene on 
benzene rings (polymer 2, which could incorporate some HCN and HC 3 N); and the polymerization 
of hydrogen cyanide and other nitriles (polymer 3). We will assume that the tholins and Titan’s 
aerosols are mainly a mixture of these types of macromolecules, to first order. 


3 Parameterization of polymer formation 

In this section, we will discuss the proposed mechanisms for polymer growth, and the necessary 
assumptions that we had to make. In each case a set of parameters will be isolated. These are 
summarized in Table II. For each pathway, initial reactions are detailed. Their rate coefficients 
( k , in cm 3 s— 1) are the first parameters. Once these reactions occured, the products are assumed 
to ultimately be incorporated in precursor particles. The parameters (a, /?, 7) characterize the 
proportions of each small molecule that is incorporated in the global macromolecule. The macro- 
molecules grow through addition of these small molecules until they reach the precursor level. 
These parameterizations allow us to calculate the mass production rates, the C/N and C/H ratios, 
and also the incorporation rates of the parent molecules (C 2 H 2 , HCN, HC 3 N and other nitriles). 

3.1 Acetylene polymer (polymer 1) 

The exact structure of the polymers of acetylene, or of a mixture of acetylene and cyanoacetylene 
(Bar-Nun ct al, 1988; Clarke and Ferris, 1997), is not well knowm. It could be chains of double 
bonds C=C, as proposed by Clarke and Ferris (1997), a family called vinylacetylenes (VA) in 
Dimitrov and Bar-Nun (1997). As a mechanism to their formation, a possibility is the attack of 
the triple bond in either acetylene or cyanoacetylene by a radical site at the end of the chain: 

Growth scheme: R n + £*£2 k B ^ak A R «+i 

where R n represents an intermediate step in the polymerization process. Its structure is shown in 
Fig. 1. This mechanism is reported as dominant for the growth of the VA family in the context of 
the model developed by Dimitrov and Bar-Nun (1997). The rates k,A and kn are not known, but 
the main parameter we will use is a = kn/kA ■ It will determine the proportions of acetylene and 
cyanoacetylene incorporated in the polymer, and the initial reactions will yield the evaluation of the 
production rates. In this parameterization, the exact process of the polymerization is not crucial. 
Once the initial products are formed, we assume that they will ultimately form precursors with a 
mass determined by the incorporated gaseous molecules. Whether other species are incorporated 
in such polymers is considered here to be second order. Clarke and Ferris (1997) described the 
incorporation of methane, ethane and carbon monoxide, and it may be possible that HCN also 
gets incorporated in those chains, but we leave this question open. 

Initial reactions: 

The initialization of the polymerization is not a simple question. It could start w r ith the radical 
C 4 H 3 or its nitrile equivalent C 4 H 2 CN: 


Initial scheme (la): 


fc ' (!) £»( 4 ) 

C2H 4- C2H2 -4 C4H2 + H — > 

C4H3, + H 

k (v) 

% 

C4H4, C4H2 + H 2 or 2 C 2 H 2 


+ (C2H2.HC3N) 

k {7) 

% 

polymer 1 (e.g. l-CgHs) 
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Table II: Proposed parameters 
Parameters 


Polymer 

1/ Poly-C 2 H 2 /HC 3 N k[ 7 ’ 9) 
a 
Ni 

2/ PAHs *4 1,2) 

Pi 

P'2 

n 2 

3/ Poly-nitriles 

7 j 


Initial reaction rate coefficients (mainly C2H2 + C4H3), in cm 3 s _1 

Ratio between HC3N and C 2 H 2 addition rates 

Mean number of molecules in a precursor macromolecule 

Initial reaction rate coefficients (C 2 H 2 + AjJ, in cm 3 s -1 

Ratio between HCN and C2H2 addition rates 

Ratio between HC3N and C 2 H 2 addition rates 

Mean number of molecules incorporated at the precursor level 

Initial reaction rate coefficient (HCNH + HCN), in cm 3 s~ l 

Ratio between nitriles and HCN addition rates 

Mean number of nitrile molecules in a precursor macromolecule 


Initial scheme (lb): 


C 3 N + C 2 H 2 
or 

k (2) u(5) 

4 hc 5 n + h 4 

C 4 H 2 CN, + H 

/.( 8 ) 
% 

c 4 h 3 cn 

C 2 H + hc 3 n 

k{ 3) 

+ (C 2 H 2 ,HC 3 N) 

uW 

% 

polymer 1 


The reactions (C 4 H 3 , C 4 H 2 CN, l-CaHs) + (C 2 H 2 , HC 3 N) — 1 polymer 1 are the first step of 

the growth scheme. 

Possible structures for C 4 H 3 and I-C 6 H 5 are shown in Fig. 1. C 3 N is the nitrile equivalent of 
C 2 H: ChC-ChN, produced by the photodissociation of HC 3 N. The reaction rate coefficients (in 
cm 3 s -1 ) k[ 7 \ k[ 4 ' and k[‘ ' have been measured (see references in Table III), but these values 
may have significant uncertainties due to different conditions between experimental works and 
Titan’s atmosphere, as discussed in photochemical models (e.g. Lebonnois et al, 2001). The rate 
coefficient k^ has been estimated in previous photochemical models, but improved data would 
be necessary. The other rate coefficients are not known. Rough estimates can be made using an 
analogy between compounds for k[ 2 \ k[ 3} , k[ 5) , A:( 8) and k[ 9 \ 

The parameter Ni (where the subscript refers to the type of polymer) is the mean number of 
molecules that are incorporated in a chain at the precursor level (by definition, the added mass 
from gas phase to the aerosol material is negligible after this yield). The relative amount of 
acetylene vs cyanoacetylene that are added depends on their respective addition rates and their 
respective number densities. The loss rate of acetylene to the formation of this polymer is 


dnc 2 n 2 

dt 


J 

kj,C 2 1 1 2 1 ^ I Ij If C ' 2 1 1 - ^fl,C 2 II 2 ) 

1=1 


( 1 ) 


where Rj ( j = 1 to J) are all the different possible radicals that could react with acetylene during 
this polymer growth, and fcj,c 2 n 2 are the reaction rate coefficients associated with these reactions. 
The loss rate is therefore related to the acetylene concentration through 

j 

^i,c 2 ii2 = y>j.C 2 ii 2 nR;, (2) 

l=i 

which we will call the ‘addition rate’ of acetylene for the formation of polymer 1. A similar term 
will be used for the different molecules and the different pathways all along this paper. The 
parameter a is defined as the ratio between cyanoacetylene and acetylene addition rates. Though 
this ratio depends in principle of the conditions in the gas under study, it may not vary much with 
these conditions if the ratio between the rate coefficient of the reaction Rj + HC 3 N and the rate 
coefficient of the reaction It- + C 2 H 2 does not depend significantly on the radical R j ; and on the 
temperature and pressure of the gas (within a reasonable range). This hypothesis is plausible, but 
certainly needs further investigation. The fractions of monomers corresponding to each molecule 
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c 4 h 3 


H-C= C-C" 

h' 


/ 

c. 


i-c 6 h 5 


+ C2H2 


H H 

/ / 

/ C ■ 

H-C= C-C O'" 

«' H 


( + C 2 H 2 or HC 3 N ) 


Rn 



H. 


(or CN\ 

/H \ /(orCN) 


(or CN) / H (or CN) 


Figure 1: Possible structures for compounds involved in polymer 1 formation. 


are: 


/l,C 2 II 2 = 
/ 1 . 11 C 3 N 


«C 2 II 2 


«C 2 II 2 + an IIC 3 N 
QJJnCsN 


( 3 ) 

( 4 ) 


nc 2 u 2 + a«nc 3 N 

where n* is the number density of the species i. The elemental ratios C/N and C/H can then 
be calculated. Using these fractions, we evaluate the mean numbers of C, H and N atoms in the 
resulting polymer chain: 

2/1.C2II2 + 3 /i : iic 3N 


(C/N)! 


(C/H)! 


/i,iic 3 n 

2/i,c 2 ii 2 + 3/i ,iic 3 n 
2/l,C 2 II 2 + /l.IICaN 


( 5 ) 


( 6 ) 


The production rate of the polymers, pi (in cm -3 s _1 ) can be evaluated as a function of the 
composition (concentrations of acetylene, cyanoacetylene and other constituents included in the 
initialisation scheme) and of the reaction rates of the initial reactions: 


n (Ta) , / (7b) , . /,( 9a) . ,(96) , 

Pi = nc i ii 3 \ K 1 ™c 2 ii 2 +«i nnc 3 N) + wc 4 h 2 cn(«i nc 2 n 2 + «i ^iic 3 n)- 
Then the mass production rate, P\ (in g cm -3 s _1 ) is: 

Pi = m. i x p l , 

where mi is the mean molecular mass of the polymer: 

mi = JVi x (/i,c 2 n 2 mc 2 ii 2 + /i,nc 3 Nmnc 3 N)- 
The loss rate of acetylene and cyanoacetylene to the polymers are also calculated: 

f dnc 2 \\ 2 


\ dt 


= Ny X /l ) C 2 II 2 X Pl, 


(7) 

(8) 
( 9 ) 

(10) 
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( 11 ) 

It should be noted that Eq. 1 and Eq. 10 are two different expressions of the same variable: Eq. 1 
is detailed from all the different reactions incorporating acetylene in the polymers, while Eq. 10 is 
obtained because (Ah/i ,C 2 n 2 ) molecules of acetylene are statistically incorporated in each chain, 
and the initial production rate of the chains is pi . 


3.2 PAH growth (polymer 2) 

If we consider a polymerization including aromatic compounds, the same formalism can be used. 
Polymerization can start through addition of acetylene on phenyl (noted A x ) , followed by an- 
other acetylene addition leading to A 2 , CioH 8 (Wang and Frenklach, 1994; Wong et al., 2000; 
Bauschlicher and Ricca, 2000). Addition would then continue and build up PAHs. The initial 
reaction considered in this pathway is: 
fc (1,2) 

A'j + C 2 H 2 A PAHs. 


HCN or HC 3 N could also be added to form heterocycles (Ricca et al., 2001). The parameters we 
will use are /?i and P 2 , where Pi (resp. p 2 ) is the ratio between HCN (resp. HC 3 N) and C 2 H 2 
addition rates. In the polymer 1, we use only two different possible monomers. Here, we have three 
different molecules that can be added to the growing macromolecule, and therefore two parameters 
(Pi, Pi) instead of one (a). Using the same formalism, the fractions of monomers included in the 
PAH corresponding to each molecule are: 


/ 2 ,C 2 II 2 - 

ftc 2 ii 2 

( 12 ) 

«c 2 ll 2 + Pl n IICN + Pin iic 3 n ’ 

/ 2 ,IICN = 

/?l».IICN 

(13) 

"'C2II2 + ‘i "lit t\ + /^ 2 » 11 G 3 N ’ 

/ 2 ,IIC 3 N = 

/? 2 «IIC 3 N 

(14) 

«C 2 II 2 + A'ftllCN +/S 2 nnc 3 N 


The initial reactions determine the production rate p 2 , then given a maximum size of grotvth (N 2 
is the mean number of molecules included in the PAH at precursor level), the mass production 
rate P 2 can be calculated: 

P -2 = m -2 X p- 2 , (15) 


with 


m-2 = Nn x (/ 2 ,C2ii2 m c 2 ii2 + h. iicnw-iicn + /2,iic 3 n?tiiic 3 n)- 


C/H is slightly over 1, depending on the size of the PAH, and C/N is 


(16) 


(C/N ) 2 


2 /2,C 2 Il2 + 3f-2,UC 3 N + /2JICN 
/2JIC3N + /2,IICN 


(17) 


The loss rates of acetylene, hydrogen cyanide and cyanoacetylene to this polymer are: 


<foc 2 n 2 

dt 


2 


{ dniicn A 

V dt ) 

( dn nCgpA 

V dt ) 


N 2 x / 2j C 2 II 2 x p 2 , 

(18) 

N 2 X /2JICN X p- 2 , 

(19) 

N 2 X / 2 ,IIC 3 N X p 2 . 

(20) 


7 



L^SJ [a] 

^N-R, 
R 1 -C 2 

z 

III 

0 

X 

+ 

— > 

^ N " R 2 

r i-cC 

1 ^C=N-H 

[b] 

^N-R, 
r i-ct 2 

+ R 3-C= N 

— > 

R-j-C^ N-R 2 

n=c 'r s 

1 M 1 [c] 

H-C= N 

+ H-C= N 

— > 

^ N-H 
H-C^ 

C=N 

[d] 

„ N-H 
or? H-C.'' 

+ H-C = N 

— 

^ N-H 

h -<c=n-h 

Initial reactions proposed 




[e] 

*N-H + 
H-C. + 

z 

III 

0 

X 

- H- 

^ N-H 

' C -C=N-H ° r H 


N-H 

N=C’ 


k <l) 

[f] H -c* N H + Rj -c= n h-c; 


N-H 

N=C' 


R. 

J 


Figure 2: Possible polymerization of nitriles (TS = Thompson and Sagan, 1989; M = Matthews, 
1992). 


3.3 HCN and nitrile polymer (polymer 3) 

The polymerization of HCN has been studied previously (Thompson and Sagan, 1989; Matthews, 
1992; Minard et al, 1998). Despite the differences in the proposed structures, it seems possible 
to reconcile them, so we will explore this possibility. Fig. 2 presents the different reactions that 
are discussed here. Thompson and Sagan (1989) describe a formation scenario initiated by the 
reaction of HCNH radical on a nitrile group. The polymer grows through subsequent N-terminus 
and C-radical interactions (reaction [b] ) . In Matthews (1992), scenarios for HCN polymerization 
proceed mainly by the way of a HCN dimer (reaction [c]). Following the discussion of Thompson 
and Sagan, it could be possible to reconcile both schemes if the first HCN dimer were formed by 
reaction [d], involving HCNH radical. This hypothesis remains to be addressed in more detail. The 
exact structure of HCN polymers remains difficult to assess, but based on these studies, we will 
make the assumption that polymerization starts with the reactions between the HCNH radical and 
nitriles (reactions [e] and [f]). In Titan’s atmosphere, the observed nitriles are HCN, then HC3N, 
C0N2, and CH3CN. We add in this study C2H3CN, which is predicted as a significant nitrile by 
our photochemical model of Titan’s atmosphere (see Section 4), and which has been detected as a 
significant nitrile in laboratory experiments (Coll et al., 1999). N 3 is the mean number of nitrile 
molecules in the precursors. Once again, the exact structure of the polymer is not crucial. The 
important factors are the initial reactions chosen, and the relative reactivity of the different nitriles. 
For each nitrile RyCN (Ry is CoH, CN, CH3 or C2H3), we define 7 y as the ratio of the addition 
rate of RyCN to the addition rate of HCN. 

The production rate of polymer 3, p 3 (cm -3 s _1 ), can be evaluated as a function of the reaction 
rate coefficients of the inital reactions and the composition. The different proportions of each nitrile 
in the polymer are linked to the relative abundance of each nitrile, and its relative reactivity: 


/3JICN 


nncN 

«IICN + Ey 7j«RjCN ’ 


( 21 ) 
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( 22 ) 


t 7j«RjCN 

/3.RjCN - — F= • 

nilCN + 7j n RjCN 

The mass production rate P 3 (g cm -3 s _1 ) is then: 


P 3 = m 3 x p 3 , 


(23) 


with 


W 3 = N 3 X (/3,IICN?R[ICN + ^ /3,RjCN«lRjCN)- 

j 


(24) 


The C/N and C/H ratios can be calculated, as well as the loss rate of the nitriles to the polymer: 


(C/N) 3 


/ 3 JICN + ") \ 

/3,IICN + f 3 , RjCN ’ 


(25) 


(C/H) s 


/3,IICN + J2j a f f 3 , RjCN 
/3,IICN + flj'/s.RjGN ’ 


(26) 


where af, 
and 


af, and a ’ 1 are the numbers of C, N and H atoms (respectively) in the RjCN molecule, 


rf'RRjCN 

dt 


- N 3 X /:;.!(,( 'X X p 3 . 
3 


(27) 


3.4 Mixing the different polymers 

The precursor level is defined as the yield where chemical growth through addition of mass from 
the gas phase becomes negligible. Chemical activity is still present (e.g. aging of the polymers, see 
Dimitrov and Bar-Nun, 2002), but the particles evolve at approximately constant mass, without 
absorbing any more gas molecules. This statement does not concern possible condensation of gas 
phase molecules on aerosol particles in the low stratosphere and troposphere (typically between 40 
and 150 km altitude), which is a problem beyond the scope of this paper. 

The total mass production rate of the haze is the sum of the different pathways, and the C/N and 
C/H mean ratios can be calculated: 

_ iVm(2/l,C2H 2 + 3 /i,hC 3 n) + ^ 2 ^ 2 ( 2 / 2 , 02 H 2 + 3 / 2 .HC 3 N + /2,HCn) + ^3P3(/3,HCN + X)j a f h, RjCN) 
NlZn/l,HC 3 X + A r 2P2(/2,HC 3 N + .D.HCn) + t\ f 3P3(/3,HCN + J2j “^/s.RjCn) 

(28) 

_ Ahpi(2/l,C 2 H 2 + 3 /i,hc 3 n) + ^ 2^2(2/2,02 H 2 + 3/2.HC3 N + /2,HCn) + A r 3P3(/3,HCN + Z)j 0^/3, RjCN) 
^1P1 (2/l,C 2 H 2 + /l,HC 3 x) + N 2 P 2 (2/o, c 2 H 2 + /2,HC 3 N + /2 ,HCn) + P3?)3(/3,HCN + Y.j a f /3,RjCN) 

( 29 ) 


3.5 Discussion of parameters 

Table II summarizes the parameters chosen for all polymers. The rate coefficients (measured or 
estimated) for the initial reactions of each scheme are given in Table III. These reaction rate 
coefficients have a significant uncertainty. The rate coefficients of HCNH on HCN (/cl 1 *) is not 
known. As a first assumption, we will estimate its value from the reaction between C 2 H 3 and 
HCN. 

Ratios between addition rates 

From a mixture of 250 torr of acetylene and 25 torr of cyanoacetylene, Clarke and Ferris (1997) 
measured a C/N ratio of 6.36, and a H/N ratio of 3.86. The polymer they report is our first 
pathway. Using Eq. 5 and Eq. 6 , and a C 2 H 2 /HC 3 N concentration ratio equal to 10, the C/N 
and C/H ratios can be fairly well reproduced for a = 6 (C/N = 6.33 and H/N = 4.33). This is 
consistent with laboratory studies showing cyanoacetylene to be two to five times more reactive 
than acetylene toward polymer formation (Clarke and Ferris, 1995, 1997). 

The studies by Ricca et al. (2001) can help constrain pi and p 3 . The addition step for HCN or 
HC 3 N in the process of ring formation show's a barrier very similar to the C 2 H 2 addition step. 
The nitrogen atom will then induce a higher barrier in the ring closure process, but this will not 
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Reaction 


Table III: Useful reaction rates 

Reaction rate coefficients (cm 3 s— 1) References 


C 2 H 2 + C 2 H -> C 4 H 2 + H 
C 2 H 2 + C 3 N -»• HC S N + H 
HC 3 N + C 2 H -x HC 5 N + H 
C 4 H 2 + H + M -x C 4 H 3 + M 

HC 5 N + H ->■ C 4 H 2 GN 

C 4 H 3 + H -> C 4 H 4 
C 4 H 3 + H -»■ C 4 H 2 + H 2 
C 4 H 3 + H -> 2 C 2 H 2 
C 4 H 3 + C 2 H 2 -> i-c 6 h 5 
C 4 H 3 + HC 3 N ( — t polymer 1) 
C 4 H 2 CN + H -t C 4 H 3 CN 
C 4 H 2 CN + C 2 H 2 ( -> polymer 1) 
C 4 H 2 CN + HC 3 N ( — t polymer 1) 
I-C 6 H 5 + C 2 H 2 ( t polymer 1) 
1-C 6 H 5 + HC 3 N ( -> polymer 1) 


feW = 8.6 x 10 -ie T i.s e 474/T 

k\ 2) = 8.6 X io -16 T 1 ' 8 e 474 / T 

fe| 3) = 8.6 X 10 - 16 T 1 - 8 e 474 / 7 ’ 

k[ A) : fcoo = 1.39 X 10 -lo e -u84 / T 

feo = 1 x 10 -28 (c) 

k[ 5) : fco„ = 1.39 X 10~ 10 e _u84 / T '! 

k 0 = 1 X 10 — 28 ? (c) 

kf a) = 8.56 x 10 - 10 e - i05 /T 

k[ 6l,) = 1.2 X 10- 11 

kf^ = 3.3 x 10- 12 

k[ 7a) = 2 x 10- 16 

k[ 7b) = 1.2 x l(r 15 

k{ S) = 8.56 X 10“ 10 e — 405 / T 7 

k[ 9a) ~ 2 X 10- 18 7 

k{ 9h) ~ 1.2 x 10- 15 ? 

2 x 10- 16 7 
1.2 x IQ" 15 7 


Opansky and Leone, 1996a 
Estimated from t) 1 * (a) 
Estimated from fc) 1 ' (a) 
Nava ef al. , 1986 
Yung ef al., 1984 
Estimated from &) 4 ' 1 (a) 


Toublanc ef al., 1995 
Yung ef al., 1984 
Yung ef al., 1984 

Wang and Frenklach, 1994 
( 7 ) 

Estimated from a X k{ 
Estimated from k\ 6a ^ (a) 
Estimated from k[ T (a) 
Estimated from ax k['^ (a) 
Estimated from fcj 7 ' (b) 
Estimated from axl)} 1 ' (b) 


Aj + C 2 H 2 AiC 2 H + H 
Aj + C 2 H 2 AiC 2 H 2 


fcj 1 ' = 6.6 x 10 17 T 1 - 58 e 1913 / T Wang and Frenklach, 1994 

k ^ = 9.8 x 10- 13 T°- 21 e- 2517 / T Wang and Frenklach, 1994 


HCNH + HCN ( ->• polymer 3) k^ = 1 x 10- 12 e- 9O °/ T ? 


Estimated from C 2 H 3 + HCN, 
Monks ef al., 1993 


HCNH + nitriles ( -> polymer 3) k^ = 7, x k { ^ parameters 

(a) If we neglect the pathway through C 4 H 2 CN against the one through C 4 H 3 (because of their relative abundance), 
then these values are secondary. 

(b) We consider that once I-C 6 H 5 is reached, all molecules will eventually yield a polymer. In this approximation, 
these rates are of secondary importance. 


(c) ko values are in cm 6 s *. 


have consequences on the / 2 factors. It may slow the growth process, but it will eventually lead 
to precursor formation. As a rough hypothesis, we consider the assumption /?i ~ /? 2 ~ 1 • 

For polymer 3, in the absence of any data, we will make the crude assumption that all nitriles have 
the same addition rate in the polymer. The parameters 7 j , ratios of the addition rate of a nitrile 
RjCN to the addition rate of HCN are therefore taken equal to 1 for all nitriles. 

Length of the polymers 

The mean length of the polymers (related to the Nj, parameters) may be controlled by the details 
of the chain reactions. To estimate the extent of the molecules’ incorporation at the precursor 
level (i.e. JV*), we use here a different approach. The usual density for aerosol particles is around 
1 g cm -3 (Toon et al, 1992; Rannou et, al, 1995). The size of the smallest particles used in 
these microphysical models for aerosol production is around 13 A. The mass of these precursors is 
therefore of the order of 9xl0 -21 g, which corresponds to ~ 200 molecules of acetylene or hydrogen 
cyanide, or ~ 400 atoms of carbon and/or nitrogen. The precursor level that is used in our model 
is taken before or around that size. Therefore, we will assume that Ni are all of the order of 10 
to 100 gas molecules (and that ~ JV 2 ~ N 3 ). Our chosen value is 20, which means that the 
precursors contains ~10 polymer macromolecules, which have incorporated each a mean value of 
20 gas phase molecules before the chemical growth from the gas phase becomes negligible. This 
is consistent with the estimation given in Dimitrov and Bar-Nun (2002), stating that a particle of 
radius (3-8) A may already consist of (12-20) macromolecules. 


10 



Table IV : Chosen parameter values 


Parameter 

Value 

a 

6 

fa, fa, 7 j 

1 

Ni, iV 2 , N 3 

20 


4 Application to the atmosphere of Titan 

4.1 Model 

In this study, we use a 1-dimensional photochemical model, updated from Toublanc et al. (1995) 
(Lebonnois and Toublanc, 1999; Lebonnois et al . , 2001). The eddy diffusion coefficient profile has 
been modified in order that the vertical profiles of the atmospheric components obtained from 
this model are similar to those obtained with our recent 2-dimensional model (Lebonnois et al., 
2001). The vertical composition of the atmosphere is sensitive to the diffusion coefficient, but 
few constraints are available: the analysis of Voyager 1/IRIS observations by Coustenis et al. 
(1989, 1995) constrain the composition in the [100-130] km region (low stratosphere), and a recent 
reanalysis of Voyager 1/UVS data by Vervack (1997) proposes vertical profiles for some compounds 
in the [500-800] km region. The 2-dimensional dynamics we introduced in the stratosphere in our 
model (Lebonnois et al., 2001) have a major impact on these profiles in this region, yielding a more 
homogeneous stratosphere. The eddy diffusion coefficient profile has been tuned here to reproduce 
this effect, and to bring most profiles in agreement with Voyager 1 observations (see Fig. 3). We 
have added to the model the production of aerosol precursors and corresponding loss of gas phase 
components, using the parameterization described previously. The chosen set of parameters is 
given in Table IV. 

Our goals here are to get a first idea of the impact of this parameterization on the photochemical 
model, to evaluate the possible production rate (its vertical profile and its column integrated value) 
of precursors in these conditions, and the composition we obtain for Titan’s aerosol precursors. 

4.2 Results 

Fig. 3 

The profiles of components of interest for this study are shown in Fig. 3: hydrocarbon molecules 
(C 2 H 2 , C 4 H 2 , benzene [A 4 ]), nitriles (HCN, HC 3 N, C 2 N 2 , CH 3 CN, C 2 H 3 CN), and radicals (C 4 H 3 , 
HCNH, phenyl [AJ). The impact on their profiles of the loss of these molecules to the aerosols 
is also indicated in this figure. As previously suggested by McKay (1996), this loss appears as a 
significant sink for the nitriles, but not for the hydrocarbons. 

The profiles obtained for the mass production rates of each polymer are shown in Fig. 4. The 
maxima are all located between 150 and 200 km altitude. This is essentially a mirror of the density 
profiles of the different radicals, the maxima being slightly lowered by the increasing densities of 
acetylene and hydrogen cyanide. The column integrated values of the mass production rates are 
indicated in Table V, at two different steps: just after the initial reactions (before addition of any 
other gas molecules), and at the precursor level. 

Polymers 1 and 3 are of similar importance, while the PAHs component (polymer 2) appears 
negligible by several orders of magnitude, due to the low mole fraction of phenyl (and therefore 
benzene) in the stratosphere, and to the addition rate of acetylene on phenyl, which is much lower 
than the other initial reactions at this temperature (in the range 150-200 K). 

The profiles of the C/N ratios are plotted in Fig. 5. This model (through Eq. 28) gives a value of 
C/N in Titan’s stratosphere of approximately 4, which is a mix of polymer 1 with (C/N) 4 ~ 18, 
and polymer 3 with (C/N) 3 ~ 1.5. The C/H ratios are always slightly above 1. 


4.3 Discussion 

Experiments 
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Figure 3: Profiles of some compounds of interest: (a) Hydrocarbons; (b,c) Nitriles; (d) Radicals. 
Solid lines: with the loss of these molecules to the haze; dashed lines: without this loss. The 
analysis of Voyager 1/IRIS spectra (Coustenis and Bezard, 1995) is indicated as boxes (or vertical 
bars for upper limits). 


Thompson et al. (1991), McDonald et al. (1994), and Coll et al. (1999) have produced tholins in 
experiments where gas phase compositions were close to Titan's stratosphere. Though the other 
experimental conditions (temperature, pressure, energy source) are not a perfect reproduction 
of Titan’s atmospheric conditions, the optical properties of tholins they produced provide good 
matches to the properties of Titan’s haze (Sagan et al, 1992). The C/N ratios they measured 
(Table I) appears to be close to the global ratio we obtain around the maximum of the production. 
The C/H ratios they report are lower than unity, which is different from our value, but this can 
be understood as a consequence of the hydrogenation of the different C-C or C-N multiple bounds 
which can occur during the polymerization process, and that is not taken into account in our 
parameterization. Another process could also be the addition of methyl groups, due to the high 
abundance of methane. 

In the Coll et al. (1999) experiment, HCN polymers (polymer 3) were reported to be likely present, 


Table V : Column integrated production rate at two different steps and the altitude of the maximum 
production rate £ max 



Initial reactions 

Precursor level (V) 

^max 


(g cm 

- 2 s- 1 ) 

(km) 

Polymer 1 

2.8xl0" 15 

2.3xl0" 14 

150-200 

Polymer 2 

3.9x 10 -19 

2.4x 10 -18 

180-220 

Polymer 3 

2.5x 10 -15 

1.8x 10 -14 

150-200 

Total 

5.3x 10 -15 

4.1 x 10 -14 

150-200 

Toon et al. (1992) 

- 

1.2x 10 -14 

200-300 

Rannou et al. (1995) 

- 

2.1 x 10 -14 

385 
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Figure 4: Profiles of the mass production rates for each polymerization pathway. 


and the C/N ratio is fairly similar to the one we get for our simulations, though slightly lower. For 
the series of experiments initiated by Khare et al. (1984), the presence of polyaromatics have been 
reported, though Thompson et al. (1991) do not indicate any production yield for benzene. It is 
worth noting that these experiments have been conducted at room temperature (300K, compared 
to approximately 170K in Titan’s stratosphere, and 100-150K in Coll et al. (1999) experiment). 
This has a large impact on the addition rate of acetylene on the phenyl radical (k^ ^ varies by 
more than three orders of magnitude between 150 and 300K), which could yield a much higher 
production rate of PAHs at room temperature, comparable to the two other kinds of polymers. 
Another possibility is that we are missing some other pathways to produce benzene (or other 
ring molecules), like the reaction proposed in Arrington et al. (1998) between 1,3-butadiene and 
metastable diacetylene. 

Microphysical models 

Two different classes of microphysical models for the haze structure and evolution have been 
proposed to reproduce several sets of observational data (Titan’s albedo, polarization of scattered 
light, high phase angle brightness). The models differ by the shape of the aerosol particles. In 
one class the particles keep a spherical shape during their growth. In the other the particles form 
aggregates of small spheres, with a fractal shape. The maxima of production and the values of 
total mass production rates needed by the Toon et al. (1992) (spherical particles) and Rannou et 
al. (1995) (fractal particles) models are indicated in Table V, for comparison with our model. 

The region of maximum mass production we obtain appears to be lower than expected by both these 
microphysical models. In our model, the position of this region is controlled by the vertical profiles 
of the radicals involved in the initial reactions (C 4 H 3 , HCNH). However, dynamical conditions 
in the stratosphere play an important role in the aerosols’ evolution and distribution, and in 
the compounds’ distribution. Winds associated with the general circulation are likely to carry 
small particles upward and may therefore result in a higher effective altitude of production than 
predicted here. This has been observed in coupled dynamics and microphysical models (Rannou 
et al submitted to Nature). 

From methane photodissociation yields, Podolak and Bar-Nun (1979) estimated the production 
rate of (b H 2 and C 2 H 4 (which is then mostly photodissociated into acetylene). Comparing this 
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Figure 5: C/N ratios for the different, pathways, and for the precursors. 

value (3xl0 -14 g cm -2 s ~ 4 ) to our production rate for polymer 1 is not straightforward since 
this polymer contains HC3N, and since the photochemical model takes into account all the other 
sinks of acetylene and ethylene. But to first, order, our value is indeed consistent with those simple 
photochemical considerations. 

The total integrated mass production rate is higher than expected by both microphysical mod- 
els, but the values we obtain at the precursor level should be taken as upper limits, since some 
intermediate molecules could be lost through condensation before they reach this level. Also, the 
influence of 2-dimensional dynamics may affect the calculated production rate, as well as the values 
expected from the microphysical models. 

To evaluate the role of the dynamical conditions in the stratosphere, it will be necessary to test 
this parameterization in a coupled model, including a general circulation model, a photochemical 
model and a microphysical model, in order to draw a more complete conclusion. 

Variations of parameters 

In order to test the sensitivity of our results against parameters variations, we varied the values 
of Ni, a, pi , p 2 , 7 j, and . First, we looked at the impact of the lengths of the chains at the 
precursor level (TV,;) on the production rates of the three kinds of polymers. Were the composition 
of the gas phase fixed, the mass production rate P t would be proportional to iV; (Eqs. 8 and 9 
for pathway 1) in the context of this model. This only means that for a given production rate 
of the first heavy molecules, the more parent molecules that are incorporated in a polymer, the 
heavier the macromolecules that constitute precursors, and the higher the total mass production 
rate of precursors. But increasing Ni can affect the vertical profiles of gas phase molecules, if this 
loss to the haze becomes a. dominant sink. Simulations were done wdth Ni = 20, 50 and 80. The 
integrated mass production rates Vi as functions of Ni are plotted in Fig. 6. Pi and P 2 vary 
significantly with N 2 and N 2 , respectively, because this sink is not dominant for the hydrocarbons, 
and therefore p\ and p 2 only slightly vary with Ni and N 2 . This is not the case for V 3 , because this 
sink is significant for the nitriles. When N 3 increases, the nitrile mole fraction in the stratosphere 
decreases, and V 3 reaches a limit. 

The parameters a, p and 7 may be dependent on the conditions in the atmosphere of Titan (or in 
the laboratory experiments). Therefore, we modified them alternatively, raising or lowering their 
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Figure 6 : Integrated mass production rate (Vi, in g cm -2 s _1 ) as a function of the parent molecules 
incorporated in the polymers at the precursor level ( N, ) . 


values by an order of magnitude. Concerning the first pathway, the value of a has an influence on 
the loss rate of HC 3 N to the haze. Variations around a = 6 does not affect CoH 2 in the gas phase, 
but has an impact on HC 3 N stratospheric abundance, since this pathway represents a significant 
loss for this compound. It also affects the (C/N)i ratio (from ~10 to ~150), as well as the mean 
molecular mass m 1 (small variations, around 5%), but the impact on the global C/N ratio and on 
the total mass production rate is small in the stratosphere (within 10%). Because polymer 2 is 
negligible compared to polymer 3, the values of Bi and /3o do not have any influence on the HCN 
and HC 3 N vertical profiles. They also have little influence on m.o and P 2 : in the case of /?i, C 2 H 2 
and HCN have similar molecular mass, and in the case of /3 2 , HC 3 N mole fraction is significantly 
lower than both other compounds. The effect of Bi and Bi variations is only visible in the (C/N )2 
ratio, but due to its negligible production rate, this does not influence the global C/N ratio. We 
tested the impact of 7 variations for HC 3 N in a similar manner. When 7 is raised for this nitrile, 
it increases its proportion in polymer 3. The molecular mass of the polymer changes accordingly 
(e.g. for HC 3 N, it raises the molecular mass), affecting the mass production rate P 3 . The (C/N )3 
is also slightly affected. The destruction rate of this nitrile increases, which affects its vertical 
profile. The overall variations of the global C/N ratio and mass production rate are small, within 
30%. 

Finally, variations on the value of the reaction rate coefficient between HCNH and HCN (C ' ) 
induce correlated variations on HCNH mole fraction profile (with only small variations on the 
nitriles’ profiles), but the production rate remains largely unchanged. This is consistent with these 
reactions being the dominant sink of the HCNH radical in this model. 


5 Conclusions 

Based on experimental data, three pathways have been proposed in order to bridge the gap be- 
tween the usual photochemical scheme and the microphysical models: polymers of acetylene and 
cyanoacetylene, polycyclic aromatic hydrocarbons (that could include some nitrogen in the rings), 
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and polymers of HCN and other nitriles. The parameterization proposed here allows us to esti- 
mate mass production rates, loss rates of gas phase molecules to the aerosols, and the C/N ratio 
obtained in the precursors for each production channel. A basic set of parameters is used to apply 
this model to Titan’s atmosphere. The following conclusions can then be derived: 

1 / We have been able to parameterize the mass production function of aerosol precursors from the 
photochemistry, given a certain number of assumptions: initial reactions leading to different types 
of polymers; relative incorporation of different kind of monomers in each polymer structure; gas 
molecules incorporation occuring until the macromolecules reach a given size (precursor level); no 
heterogenous nor ion chemistry considered. 

2/ Though the comparison of the resulting C/N ratio and main products to the experiments 
are promising, the values determined for the mass production rate and the altitude of the peak 
production are only marginally in agreement with the values that 1 -dimensional microphysical 
models require. 

3/ Variations of the different parameters around their chosen values do not significantly affect the 
results, except for the length of the polymers at the precursor level. 

4/ This approach can help bridge the gap between gas phase molecules and the macromolecules 
that are building the aerosols in a simple way, without getting into the details of the polymers 
structures. For this parameterization to be effective, it needs a qualitative understanding of the 
initial reactions, as well as a quantitative assessment of the different parameters. Efforts should 
continue in this direction, in particular through more detailed analysis of experimental data. The 
simultaneous characterization of the tholins (types of bounds, aliphatics vs PAHs, HCN polymers, 
C/N and C/H ratios) and measurement of the gas phase composition (mole fractions of as many 
compounds as possible) are needed to help constrain the production mechanism of the tholins. 
It would be useful to have such experimental data in different initial conditions, especially with 
different initial mixtures: apart from the usual N 2 /CH 4 , it could also be C 2 H 2 /HCN (14/1), for 
example. 

5/ The first results we obtain here should also lie tested with a coupled model (dynamics, photo- 
chemistry and haze microphysics), especially since dynamics can affect significantly the required 
production conditions. This will help build a complete model of Titan’s atmospheric system, as a 
powerful tool for the exploitation of Cassini/Huygens data. 
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